Nuclear modification of heavy quark fragmentation function and J /ip production in 

ultrarelativistic heavy ion colhsions 



Bin Zhang, Bao-An Li, Andrew T. Sustich, and Charles Teal 
Department of Chemistry and Physics, Arkansas State University, 
P.O. Box 419, State University, Arkansas 72467-0419, USA 
(May 22, 2002) 

In ultrarelativistic heavy ion collisions, charm quark fragmentation is modified due to the melting 
of strings inside the high density partonic matter. D mesons produced in hadronization can further 
produce J/tp particles. Using a multiphase transport model, we investigate the effect on the rapidity 
distribution of the produced J/ip particles with two difi'erent charm quark fragmentation functions. 
It is shown that the J/tp rapidity distribution is sensitive to the nuclear modification of charm 
quark fragmentation and thus is a good indicator of the onset of string melting and the production 
of deconfined partonic matter. 
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One interesting topic in strong interaction physics that has attracted much attention is the study of the phase 
transition between hot and dense hadronic matter and the Quark-Gluon Plasma (QGP) The QGP is believed to 
have existed in the early universe. It may also exist in the cores of neutron stars. Knowledge about the QGP phase 
transition is critical to the understanding of strong interaction physics and several other processes in cosmology and 
astrophysics. Several accelerators have been built and many experimental observables have been proposed for the 
search of the QGP J/ip production is one of the most promising probes for the formation of the QGP ||^-|2^. In 

this Letter, we demonstrate that due to the melting of strings inside a QGP, the heavy quark fragmentation function 
is modified. This modification is reflected clearly in the J/ip rapidity distribution. The measurement of the latter 
thus can be used to study the onset of string melting and the production of the QGP. 

In relativistic nuclear collisions, the charm quark distributions can be calculated from perturbative quantum chro- 
modynamics (pQCD). After the parton evolution, these charm quarks are converted into D mesons. The D meson 
distributions are determined by the charm quark distributions and the fragmentation function. In hadron induced 



reactions, charm quark fragmentation can be described by the delta function fragmentation function 23 1, 

D^iz)=Siz-l). (1) 
Charm quark fragmentation in e'^e~ annihilation can be well described by the Peterson fragmentation function [p4[ , 

- .(1-1/.-" /(I -.))' - 

In the above equations, z is the fractional momentum of the D meson relative to the charm quark. is a normalization 
factor, and Ec ^ {^q/i^cf' ■ An explanation of the observed harder D meson distribution in hadron induced interactions 
is provided by the Lund fragmentation model. In the Lund model based event generator, PYTHIA |26|, the 
charm quark is always at the end of a string and can be pulled by the faster valence particles to acquire a harder 
distribution. 

In nuclear collisions, two scenarios can happen. In scenario 1, the produced charm (c) and anticharm (c) quarks are 
pulled by the strings of the receding nuclei. These strings break up after being stretched. A c quark then combines 
with a light antiquark (g) to form a D meson. A similar process happens for the formation of the D meson (Figure^). 
In scenario 2, the c and c quarks are produced together with many gluons {g). The interactions between the charm 
quarks and the receding nuclei are screened, and the strings between them melt |2^,^^. At a later time, gluons 
fragment into quark-antiquark pairs, while other quark-antiquark pairs are also produced by the wounded nucleons. 
A D meson can be formed by the combination of a charm quark and a light antiquark (Figure ||). 

These two scenarios can lead to very different rapidity distributions for the produced D and D mesons. In scenario 
1, the charm quarks are pulled by the faster receding hadrons. The rapidity distribution will be harder and wider 
compared with that of scenario 2 in which the strings that pull the charm quarks are absent. In the second scenario, 
the charm quarks are also slowed down by other quarks, and the rapidity distribution will be softer and narrower. 
Charm quark fragmentation in scenario 1 can be modeled by the delta function fragmentation while that in scenario 
2 can be modeled by the Peterson fragmentation function. 

In ultrarelativistic nuclear collisions, the initial stage is dominated by minijet gluons. In particular, for Au+Au 
collisions at 2QQA GeV, the pQCD calculated gluon rapidity density is about 300 at midrapidity p^. For comparison. 
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the charm and anticharm rapidity density is about 3 ||29[]. In other words, these charm and anticharm quarks are 
immersed in a "sea" of hot gluons. UnHke for nuclear colhsions at lower energies, or hadron-hadron collisions, the large 
amount of gluons modify the interactions between the charm quarks and the valence quarks. The strings between 
them melt when the gluon density is high. Several calculations have demonstrated the possibility of string melting 
at the Relativistic Heavy Ion Collider (RHIC) energies j2^,^. The melting of strings then leads to the absence of 
the pulling force that is responsible for the harder D meson distribution observed in hadron induced reactions. This 
modification of the fragmentation function of charm quarks therefore reflects the production mechanism of the parton 
matter in ultrarelativistic nuclear collisions. 

We propose to use the J /if) rapidity distribution as a measure of this modification. J /ip production is affected by 
several factors. Due to color screening inside the QGP, J/ij] particles are dissociated above a certain temperature. 
A J /ip can be produced by the collision of a c quark and a c quark when the temperature is below the dissociation 
temperature. In the parton phase, a J /ip can be destroyed by its collision with a minijet gluon. Similar processes 
can occur in the hadron phase. A J /ip can be produced by the collision of a D meson and a D meson. It can be 
dissociated by the collision with a meson that is made of light quarks. The above J /ij} production from final state 
hadronic interactions relates the J /ip rapidity distribution closely to the D and D meson distributions. 

The modification of charm quark fragmentation and its reflection in the J /ip rapidity distribution are studied by 
using a multiphase transport model [^ - |32| . This model is a hybrid model that is made up of several Monte Carlo 
models. The heavy ion jet interaction generator (HIJING) [^,0 is used to simulate the passing of two heavy nuclei 
and the production of the initial parton system. Zhang's parton cascade (ZPC) |35| is then used to study the evolution 
of this parton system. The conversion of the parton matter into hadron matter when no more parton interactions 
are available is done using the Lund fragmentation |^ that is part of the PYTHIA/JETSET Monte Carlo package 
pGf . A relativistic transport (ART) model will follow the hadron evolution until freeze-out. The multiphase 

model incorporates the major channels in different stages. It has been used to study particle production at RHIC 
energies and can successfully describe the global feature of the newly available RHIC data. 

For J /ip production, c + c J /ip + is included in the parton stage and D + D ^ J /ip + M is included in the 
hadron stage. Here M is a meson that is composed of one light quark [q) and one light antiquark (q). In relativistic 
nuclear collisions, charm quarks are rare particles. To effectively simulate the evolution of charm particles and J/?A's, 
the perturbative method |^,^ is used. In this approach, charm quarks from several independent events are allowed 
to evolve with a common hadron environment. Since charm quarks are rare particles, their effects on the common 
hadronic environment are neglected. 

We will compare the J/?/; rapidity distributions in relativistic nuclear collisions for the two fragmentation scenarios. 
The multiphase transport model is used to simulate central Au+Au collisions at RHIC at a center of mass energy, 
— 200A GeV. The PYTHIA model is used to generate the initial charm and anticharm quarks. The charm 
production cross section per nucleon-nucleon collision is chosen to be 350 /zb. The mass of the charm quark is 
chosen to be 1.50 GeV, and the mass of the D mesons is set to be 2.01 GeV. A 3 mb cross section is used for the 
elastic collisions between gluons and charm quarks and for the J /ip dissociation by gluons and light mesons. The 
cross sections of the inverse processes are calculated by using the detailed balance relation. These cross sections are 
consistent with previous pQCD and hadron model calculations p3| , ^ -^ . 

The rapidity distributions of the J/ip particle, together with those for the charm quarks, D and D mesons are 
calculated. Figure |^ shows the results for the case in which the delta function fragmentation is used to turn the 
charm quarks into D mesons. It is seen that initial charm quark rapidity distribution has a plateau with a value of 
about 3.6 around y = Q. After the parton evolution, the final charm quark rapidity distribution is almost the same 
as the initial distribution. This distribution leads to a J /ip rapidity distribution that also has a plateau with a value 
of about 0.7/300 around y = 0. The rapidity distribution for the D and D mesons produced from hadronization is 
similar to the charm and anticharm quark distribution. The value around y = is about 4.3 and remains almost 
unchanged during the hadron evolution. J /ip particles are produced further in the hadron stage and the final rapidity 
distribution follows that of the D and D mesons. The final value of J /ip particle around y = is about 1.9/300. 

The distributions for the Peterson fragmentation case are shown in Figure ^. The initial, final charm quark 
distributions and the J /ip distribution after the parton stage are the same as those for the previous case since the 
parton stage has no information about the hadronization. As expected, the D and D rapidity distribution becomes 
narrower by about 0.6 units of rapidity and has a high peak with a value about 6.4 aX y — Q. Notice that the total 
number of charm quarks does not change during hadronization. Hence, the area under the rapidity distribution of 
the charm and anticharm quarks equals that of the D and D mesons. However, the production of J /ip particles is 
now starting from a D and D rapidity distribution that is different from that of the delta function fragmentation 
case. This leads to a final J /ip distribution that is also narrower and has a peak in the middle with a value of about 
5.6/300 around y = 0. 
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In conclusion, the production of dcconfined matter modifies the fragmentation of the charm and anticharm quarks. 
In nuclear collisions at lower energies, the string dynamics is important, and the charm fragmentation is similar 
to that in the hadron collisions. In nuclear collisions at ultrarelativistic energies, the minijet dynamics becomes 
important, and the charm fragmentation is more like the independent fragmentation in electron-positron collisions. 
This modification leads to a rapidity distribution of J /ij) distinctively different from what is normally expected in 
hadronic reactions. The J/ij) rapidity distribution measurement, in conjunction with information about charmed 
mesons, can offer valuable insight about the melting of strings and the production of the QGP in relativistic nuclear 
collisions at RHIC. 
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FIG. 1. Time evolution of the collision in fragmentation scenario 1. (a) Two heavy nuclei at relativistic speeds approach 
each other, (b) After the collision, a charm quark and an anticharm quark are produced. They are attached to the receding 
nuclei by strings, (c) At a later time, quark-antiquark pairs are produced from the strings, (d) The charm quark and a light 
antiquark combine to form a D meson, while the anticharm quark and a light quark combine to form a D meson. 




FIG. 2. Time evolution of the collision in fragmentation scenario 2. (a) Two relativistic nuclei approach each other, (b) 

A charm quark and an anticharm quark arc produced together with a large number of gluons. (c) At a later time, light 
quark-antiquark pairs are produced, (d) The charm quark and a light antiquark combine to form a D meson, while the 
anticharm quark and a light quark combine to form a D meson 
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FIG. 3. Rapidity distributions of J/?/" particles, charm and anticharm quarks, D and D mesons for the deha function 
fragmentation case. The dotted hne J/^l) distribution is 300 times the distribution after the parton stage. The soUd line J /^j) 
distribution is 300 times the distribution after the hadron evolution. The dotted line charm and anticharm quark distribution is 
before the parton evolution and the solid line is after the parton evolution. The dotted line for the D and D meson distribution 
is before the hadron evolution and the solid line is after the hadron evolution. 




FIG. 4. Rapidity distributions of J/t/; particles, charm and anticharm quarks, D and D mesons for the Peterson fragmen- 
tation function case. The symbols are defined the same as those in Fig. |^. 
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